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AFIT/GAE/ENY/04-M04 


Abstract 


The  purpose  of  this  study  was  to  determine  the  feasibility  of  launching  a  vehicle 
based  on  the  Boeing  AirEaunch  System  in  a  coplanar,  direct  to  rendezvous  trajectory 
with  gravity-turn.  The  focus  of  the  research  was  to  model  the  launch  trajectory  and 
determine  the  ability  to  reach  different  coplanar  orbits.  The  launch  trajectory  was 
modeled  using  two-dimensional  equations  of  motion  and  a  boundary  value  problem  was 
posed  and  solved  for  the  gravity-turn  trajectory.  Trajectories  were  then  created  in  an 
attempt  to  reach  different  altitudes  through  coasting  and  transfer  orbits.  Einally  a  specific 
orbital  altitude  was  chosen  and  the  trajectories  were  analyzed  to  find  the  most  efficient 
route  to  the  target  orbit  for  fuel  and  time. 
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COPLANAR  AIR  LAUNCH  WITH  GRAVITY-TURN 


LAUNCH  TRAJECTORIES 

I.  INTRODUCTION 


1. 1  Motivation 

The  goal  for  this  research  was  to  investigate  launching  a  vehicle  direct  to  rendezvous 
with  an  orbiting  satellite  using  a  gravity-turn  trajectory.  Current  launch  platforms  use 
launch  windows  for  mission  planning.  This  practice  launches  the  vehicle  into  the  desired 
orbital  plane  as  it  passes  overhead.  Once  on  orbit,  the  vehicle  then  has  to  play  catch  up 
with  the  target.  In  some  cases  the  catch  up  period  can  take  up  to  three  days.  The  launch 
windows  themselves  are  very  restrictive  and  sometimes  offer  limited  time  for  a 
successful  launch.  The  trade  off  for  this  restrictive  schedule  is  a  savings  in  fuel  and 
energy  and  provides  maximum  payload  to  orbit. 

Currently  the  Air  Force  is  investigating  concepts  to  provide  rapid,  responsive  access 
to  space.  Several  new  designs  attempt  to  achieve  these  concepts  through  different  means. 
One  such  possible  design  is  air  launch.  An  air  launch  system  would  be  capable  of  flying 
to  a  desired  point  for  launch,  thus  giving  the  ability  to  reach  launch  points  at  the  user’s 
discretion.  The  desired  capability  of  rapid  access  to  space  with  the  ability  to  launch 
direct  to  rendezvous  would  offer  a  range  of  options  such  as  re-supply,  emergency  repair. 


constellation  regeneration,  satellite  protection  and  many  other  opportunities  desired  by 
the  modern  military. 

Gravity-turn  trajectories  are  also  an  item  of  current  interest.  Launeh  vehicles  must 
maintain  a  zero  angle  of  attack  during  launch  through  the  atmosphere  due  to  structural 
strength.  Even  a  small  angle  of  attaek  ean  mean  structural  failure  for  the  vehicle.  In  a 
gravity-turn  trajeetory,  the  vehicle  takes  advantage  of  the  force  of  Earth’s  gravity  in  order 
to  rotate  from  vertical  to  a  horizontal  flight  orientation  tangential  to  its  orbit.  This  allows 
the  vehicle  to  eonserve  fuel  and  the  mass  of  extra  engines.  In  a  gravity-turn,  roll  and 
angle  of  attack  are  maintained  at  zero,  so  that  no  lift  is  generated.  Space  launeh  vehicles 
are  made  to  be  very  strong  along  their  longitudinal  axis,  however  are  very  weak  along  the 
lateral  axis.  If  lift  is  generated  the  vehicle  will  more  than  likely  disintegrate.  To  avoid 
this,  the  vehicle’s  computer  will  eompensate  to  keep  the  angle  of  attack  and  roll  zeroed 
out  while  letting  the  earth’s  gravity-turn  the  vehiele.  The  vehiele  is  given  a  very  small 
nudge  from  the  vertical  to  begin  the  proeess.  During  this  time  a  small  amount  of  lift  will 
be  created,  so  the  proeess  is  begun  shortly  after  launeh  when  the  vehicle’s  speed  is  very 
slow.  It  cannot  be  done  from  the  initial  launch  position  as  the  vehiele  does  not  have 
enough  momentum  and  will  simply  fall  over.  Eike  generating  lift  at  speed,  this  is  a  very 
bad  situation  for  a  rocket  to  find  itself. 


1.2  Overview 

The  rest  of  this  chapter  is  devoted  to  discussing  the  Boeing  AirEaunch  vehiele 
used  as  the  basis  for  the  researeh.  Chapter  2  deseribes  setting  up  the  trajectory  problem 
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and  introduces  the  derivation  of  the  equations  of  motion  and  Hohmann  transfer.  The 
computer  algorithms  used  in  the  researeh  are  discussed  in  Chapter  3.  Results  of  the 
gravity -turn  trajectory  analysis  are  presented  in  Chapter  4.  Chapter  5  then  finishes  with 
the  conclusions  and  recommendations. 

1.3  Vehicle  Background 

One  of  the  most  desired  eapabilities  of  eurrent  space  launch  is  the  ability  to 
launch  on  demand  into  a  first  pass  orbit.  Current  systems  sueh  as  Expendable  Launeh 
Vehicles  (ELVs)  and  the  Shuttle  Eaunch  System  (SES)  are  based  on  timetables, 
schedules,  and  launch  windows.  These  systems  require  months  and  possibly  years  of 
planning  and  executing  for  a  launch.  The  only  launch  on  demand  systems  eurrently 
operating  is  that  of  Orbital  Seiences  of  Dulles,  Va.  In  Orbital’s  system  a  Pegasus  launeh 
vehiele  is  carried  aloft  on  a  Eoekheed  E-101 1  and  flown  to  a  speeific  point  along  the 
launch  corridor.  The  Pegasus  can  lift  small  payloads  of  1,000  lb  into  orbit.  Boeing’s 
proposed  AirEaunch  System  is  an  attempt  to  support  payloads  on  the  order  of  15,000  lb. 
The  system  will  be  carried  on  the  back  of  a  Boeing  747-400E  freighter  and  have  the 
ability  to  operate  from  any  10,000  to  12,000  foot  runway.  (Wilson,  2001:43-46) 

The  Boeing  747-400E  would  then  carry  the  launch  vehicle  to  a  speeifie  launeh 
altitude  and  position.  With  an  in-flight  refueling  capability,  practically  any  launch  point 
can  be  reaehed.  AirEaunch’s  wings  and  tail  assembly  will  provide  the  vehiele  with  glide 
ability  so  that  a  safe  distance  ean  be  aehieved  between  separation  from  the  aircraft  and 
engine  ignition.  Approximately  five  seconds  after  ignition  the  vehicle  would  jettison  its 
wings  and  tail  for  the  launch.  It  should  be  noted  that  the  discussion  on  operations  from 
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Boeing  does  not  give  a  great  amount  of  detail  on  the  length  of  the  pull  up,  or  the  lateral 
aeeeleration  tolerable  by  the  vehicle.  For  the  purposes  of  this  research,  it  is  assumed  that 
the  vehicle  can  reach  a  vertical  attitude  safely  at  20,000  feet  and  a  velocity  of  300  mph. 
("Phantom  Works”,  2003) 

Boeing’s  AirLaunch  System  is  the  basis  of  the  vehicle  used  in  this  research.  The 
launch  vehicle  of  AirLaunch  will  consist  of  3  stages  plus  payload.  The  first  two  stages 
are  Thiokol  Castor  120’s  with  the  third  stage  made  up  of  a  Thiokol  Star-92.  Castor  120’s 
are  off  the  shelf  solid  rocket  motors  with  specifications  shown  in  Table  1. 


Table  1:  Castor  120  Specifications  (“Castor  120”,  2003) 


Length 

347  in 

881.38  cm 

Diameter 

93  in 

236.22  cm 

Propellant  Weight 

107,137  Ibm 

48,596.526  kg 

Total  Weight 

116,159  Ibm 

52,688.836  kg 

Average  Thrust 

370,990  Ibf 

1,650,245.737  N 

Specific  Impulse 

280  sec 

280  sec 

Bum  Time 

82  sec 

82  sec 

Thiokol  Star-92 ’s  solid  rockets  are  not  currently  released,  and  no  information  was 
available  from  Thiokol  or  Boeing.  Therefore,  specifications  were  assumed  from  known 
solid  rocket  motor  performance  and  limitations.  The  ratio  of  structure  mass  to  payload 
mass  was  also  assumed  to  be  similar  to  those  of  the  Castor  120’s.  A  mass  flow  half  that 
of  the  Castor  120  was  used  to  produce  the  burn  times  and  average  thrust.  For  the  mass  of 
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the  third  stage,  the  given  total  vehiele  weight  was  used  minus  the  known  payload  mass 
and  the  masses  of  the  first  two  stages.  The  specifieations  used  for  the  Star-92  stage  is 


given  in  Table  2. 

Table  2:  Star-92  Specifications 


Propellant  Weight 

25,177.853  kg 

Total  Weight 

27,298.096  kg 

Average  Thrust 

813,932.545  N 

Specific  Impulse 

280  sec 

Bum  Time 

84  sec 

The  major  motivation  to  use  the  Boeing  AirLauneh  System  as  the  basis  for  this 
research  was  for  the  capabilities  of  air  launch.  The  focus  of  the  research  was  to  look  into 
the  feasibility  of  a  gravity-turn  launch  trajectory  into  a  direct  to  rendezvous  orbit  for  the 
Air  Force’s  Space  Maneuvering  Vehicle  (SMV).  To  achieve  this  goal  with  launch  on 
demand,  a  movable  launch  pad  or  air  launch  was  preferred.  Boeing’s  AirLauneh  is  the 
closest  feasible  vehicle  in  achieving  the  ability  to  launch  on  demand  in  the  payload  class 
of  the  SMV  of  7,500  lb  (3,000  kg). 


5 


II.  MODELING  THE  TRAJECTORY 


2.1  Literature  Search 

A  search  of  current  literature  produced  articles  on  the  Boeing  AirLaunch  System 
(ALS).  This  information  included  stage  types,  total  weight  of  the  vehicle,  specifications 
on  the  first  and  second  stage  engines,  and  some  operational  information.  The  ALS  is 
piggy  backed  aloft  on  a  Boeing  747.  The  operations  of  this  carrier  craft  were  discussed, 
allowing  types  of  air  fields  capable  of  deploying  from,  and  an  in  flight  refueling 
capability  giving  the  carrier  aircraft  a  virtual  unlimited  range.  The  vehicle  can  then  be 
flown  to  a  launch  point  as  necessary.  However,  no  trajectory  profiles  were  discussed  in 
the  research. 


2.2  Reference  Frames 

The  frame  decided  upon  for  this  research  was  simplified  to  the  downrange  distance 
measured  in  the  x-direction  and  altitude  measured  in  the  y-direction.  As  seen  in  Figure  1, 
the  H  axis  is  oriented  to  point  away  from  the  center  of  the  Earth,  and  the  frame  rotates 
such  that  H-axis  remains  pointed  away  from  the  center  of  the  Earth.  The  X-axis  is 
defined  as  shown  perpendicular  to  the  H-axis  in  the  direction  of  movement. 


6 


This  rotation  is  dependent  on  the  horizontal  veloeity  of  the  frame,  sueh  that; 


CO 


X 


iK+H) 


A 

-h. 


(2-1) 


Where  is  the  radius  of  the  Earth,  and  H  is  the  altitude  of  the  vehiele. 


X  is  the 


veloeity  in  the  X-axis,  and  would  be  represented  by  the  plane  eoming  out  of  the  paper. 

The  rotation  is  negative  as  it  is  opposite  of  the  standard  right-handed  notation  (Wiesel, 
1989:19). 
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The  centripetal  acceleration  can  then  be  found  as: 


o*'  x((y*'x  r)  = 


(K  +  H) 


hi 


(2-2) 


where  r  =  (Rg  +  H)  /22  and  /?2  is  in  H-axis  as  drawn  in  the  Figures  1  and  2.  However,  it  is 
conducive  in  this  situation  to  use  the  body  frame  in  the  derivations  of  the  equations  of 
motion,  so  the  flight  path  of  the  vehicle  can  be  taken  into  account  (Wiesel,  1989:217). 

As  shown  in  Figure  2  the  flight  path  angle,  y,  is  the  angle  between  the  local  horizon  and 
the  velocity  vector. 


Figure  2:  Reference  frame  with  flight  path  angle 


With  the  frames  defined,  we  are  ready  for  the  equations  of  motion. 
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2.3  State  Variables  and  Equations  of  Motion 

The  development  of  the  equations  of  motion  here  is  a  modified  form  of  those 
presented  in  (Wiesel,  1989:  216-219).  First,  the  state  spaee  vector  must  be  defined.  The 
five  elements  used  for  this  research  are: 

m  -  mass  of  the  vehicle 

H  -  altitude  of  the  vehicle 

X  -  distance  downrange  of  the  vehicle 

V  -  velocity  of  the  vehicle 

Y  -  flight  path  angle  of  the  vehicle 

With  the  state  vector  defined,  the  equations  of  motion  can  be  developed  further. 
Using  Figure  3  as  a  reference,  we  can  make  some  observations.  First,  the  vehicle  stays 

on  the  H  axis,  therefore  the  vertical  acceleration  is  H  and  the  downrange  acceleration  is 
represented  as  X  .  By  geometry,  the  following  equations  can  be  produced: 


dX 

—  =  V  cosy 

(2-3) 

dt 

dH  .  .  . 

- =  F  sm  y 

dt 

(2-4) 
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X- 

down  range 

Figure  3:  Diagram  of  Forces  Acting  on  the  Vehicle 


As  mentioned  in  the  previous  seetion,  it  is  more  eondueive  to  use  the  body  frame 
instead  of  resolving  into  the  vertieal  and  horizontal  frames.  By  doing  so,  the 
aeeelerations  of  the  vehiele  ean  be  summed  up  into  aeeeleration  along  and  transverse  to 
the  vehiele’s  axis.  Considering  that  V  is  the  veloeity  of  the  vehiele,  and  y  is  the  flight 

dV 

path  angle  measured  eountereloekwise  from  the  horizon,  then - is  the  aeeeleration 

dt 


produeed  along  the  axis  of  the  vehiele,  and  V 


dy 

dt 


is  the  aeeeleration  transverse  to  the 


vehiele’s  axis. 


Now,  using  Newton’s  Seeond  Law: 


=  ma 


(2-5) 
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We  can  then  proceed  to  sum  the  forces  of  the  vehicle.  Again  referencing  Figure  3,  we 
find  the  forces  acting  on  the  vehicle  are  Thrust  (T),  Drag  (D),  gravity  (mg),  and 
centripetal  acceleration  (2-2).  Rotating  the  gravity  and  centripetal  accelerations  by  the 
flight  path  angle  of  the  vehicle  and  summing  the  forces  along  the  axis  of  the  vehicle  then 
gives: 


-ma  -  T  - D  -mgsiny  + 


mX^ 

iK+H) 


sin/ 


simplifying  to: 


(2-6) 


dV 
m  — 
dt 


T  -D- {mg 


mX^ 

{Re  +  H) 


)sin/ 


(2-7) 


Then  similarly,  summing  the  forces  transverse  to  the  axis  of  the  vehicle  produces: 


=  ma  =  -mg  cos  /  + 


mX^ 

(Re  +  H) 


cos/ 


(2-8) 


which  then  simplifies  to: 


mV 


dy 

dt 


-{mg 


mX^ 

(Re+H) 


)cos/ 


(2-9) 
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To  find  the  mass  fiow  rate  of  the  vehiele  the  following  equation  was  required: 


m  =  -- 


g  I 

o  sp 


(2-10) 


where  T  is  thrust,  Isp  is  the  specifie  impulse  of  the  engine,  and  go  is  gravitational 
aoeeleration  at  sea  level.  For  the  third  stage,  equation  (2-10)  holds  true,  however  it  does 
not  do  so  for  the  first  and  seeond  stages.  As  the  exaet  speeifieations  are  available  for  this 
solid  rocket  motor,  a  double  check  reveals  that  with  a  burn  time  of  82  seconds,  the  mass 
fiow  provided  by  (2-10)  results  in  a  total  fuel  bum  that  is  greater  than  the  actual  available 
fuel.  Therefore,  in  the  first  two  stages,  the  mass  fiow  is  determined  by: 


total  mass  of  fuel 
burn  time 


(2-11) 


While  it  is  most  likely  tme  that  there  is  still  some  fuel  left  before  the  stage  is  jettisoned 


and  the  fuel  is  probably  not  burned  at  a  constant  rate,  this  provides  a  reasonable  model. 
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Therefore,  to  form  the  sealar  equations  of  motion  for  the  system,  we  gather  the  equations: 


dX 

dt 


V  cosy 


dH 

dt 


V  siny 


dV 

m - 

dt 


T  -D-  {mg 


mX^ 

{R.+H) 


)siny 


(2-12) 


mV^  = 
dt 


-{mg 


mX^ 


iR+H) 


)cosy 


T 

m  - - 

S  I 
o  sp 

In  general,  the  equations  of  motion  for  non-linear  time-dependent  systems  is 
written: 


X  =  f{x,u,t)  (2-13) 

Given  x  represents  the  state  variables,  u  is  representative  of  the  control  variables  and  t 
represents  time  (Sears,  1997:14-15).  This  form  is  used  by  MATLAB,  which  is  discussed 
in  the  following  chapter. 
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2.3  Hohmann  Transfer 


As  will  be  discussed  in  Chapter  4,  several  options  are  available  for  achieving  orbit. 
One  possible  augmentation  to  the  launch  is  the  addition  of  a  transfer  orbit  to  insert  the 
vehicle  into  the  proper  orbit.  The  Hohmann  Transfer  is  the  most  efficient  use  of  available 
fuel  with  the  disadvantage  of  being  the  longest  transfer  available  in  terms  of  time. 

As  derived  and  shown  in  (Wiesel,  p74-75),  there  are  two  changes  in  velocities  or 
AV’s  required.  These  two  velocity  changes  are  represented  by  equations  2-14  and  2-15. 


AV, 


2ju  Iju 


(2-14) 


(2-15) 


where 

AVi=  change  in  velocity  at  first  maneuver 
AV2=  change  in  velocity  at  second  maneuver 
q  =  gravitational  parameter  of  earth 
ai=  radius  of  the  first  maneuver 
a2=  radius  of  the  second  maneuver 
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The  gravitational  parameter  of  earth  (p)  is  defined  as  3.98601  x  10^  kmVs^  (Wiesel, 
1989:323).  A  physical  relationship  between  the  variables  in  equations  2-14  and  2-15  can 
best  be  seen  in  Figure  4  below: 


Figure  4:  Hohmann  Transfer  between  two  orbits  (Wiesel,  1989:74) 
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Also  of  interest  is  the  time  of  flight  of  the  transfer  orbit.  A  Hohmann  transfer  is  basically 
half  of  an  elliptical  orbit.  So  as  seen  in  (Wiesel,  1989:75)  the  time  between  maneuvers  is 
given  as: 


(2-16) 


where 


a,  +  a, 

a  =  — - ^ 

2 

p  =3.98601  X  10^  kmV 


(2-17) 


The  time  between  maneuvers  is  especially  of  interest  for  mission  planning. 
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III.  ALGORITHMS 


3.1  Introduction 

In  previous  chapters  the  vehicle  and  equations  of  motion  were  discussed  and 
developed.  In  this  chapter  we  discuss  the  algorithms  used  in  the  computer  programs  to 
model  the  launch  trajectories  of  the  vehicle.  Examples  of  the  m-files  can  be  found  in 
Appendix  B.  Before  the  algorithms  can  be  truly  discussed,  the  assumptions  for  this 
problem  need  to  be  stated. 

1 .  Drag  is  being  neglected,  i.e.  D=0. 

2.  Vehicle  specifications  are  as  stated  previously. 

3.  Third  stage  assumptions  are  accurate. 

4.  The  vehicle  has  pitched  to  the  full  vertical  position  without  expenditure  of  fuel 
at  launch. 

5.  The  vehicle  begins  with  a  300  mile  per  hour  velocity  at  20,000  feet. 

6.  The  derivations  and  assumptions  stated  in  the  preceding  chapters  are  assumed 
correct. 

The  programs  were  written  in  the  following  units: 

Mass  =  kilograms  (kg) 

Time  =  seconds  (s) 

Distance  =  meters  (m) 

Velocity  =  meters/second  (m/s) 
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With  the  above  assumptions,  computer  programs  were  written  to  solve  the 
boundary  value  problem  for  the  gravity-turn  equations  of  motion.  These  programs  were 
written  in  MATLAB  6.5  Student  Edition  and  are  commonly  referred  to  as  m-files  in  this 
research.  The  entire  trajectory  profile  boundary  value  problem  was  broken  into  three 
parts.  Each  section  represents  a  stage,  as  there  is  a  discontinuity  at  the  separation  points. 
Each  program  consists  of  two  m-files.  The  first  file,  INITIALCOND  provides  initial 
conditions  and  the  initial  state  vector.  The  MATLAB  command  ODE45  then  calls  upon 
the  second  m-file,  LAUNCHEOMS  in  this  case,  that  contains  the  equations  of  motion. 
The  routine  continues  for  a  set  time  period  providing  an  output  file  of  the  state  vector  at 
each  time  step.  The  following  stage’s  m-file  then  reads  the  last  previous  state  vector  and 
runs  an  almost  identical  algorithm,  the  major  changes  being  the  initial  conditions  and  in 
the  case  of  the  third  stage,  the  thrust.  The  end  result  is  an  output  file  with  final 
conditions  in  mass,  altitude,  downrange  distance,  velocity,  and  flight  path  angle  of  the 
vehicle.  The  output  file  also  contains  the  state  vector  at  each  of  the  iterations,  providing 
data  for  manipulation  to  be  presented  in  the  following  chapter.  Other  programs  were 
used,  but  will  be  discussed  later  in  the  chapter.  The  purpose  of  these  algorithms 
ultimately  is  to  find  the  altitude,  velocity,  and  initial  flight  path  angle  that  will  provide  a 
zero  flight  path  angle  at  burn  out. 


3.2  Initial  Launch  Trajectory 

An  initial  trajectory  was  produced  first.  Several  runs  were  made,  and  typically 
the  first  trajectory  of  each  run  was  made  with  a  flight  path  angle  of  89.5°.  This  provided 
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a  reference  trajectory  for  each  run.  Obviously  an  initial  flight  path  angle  of  90”  would 
not  be  useful  as  there  would  not  be  a  gravity-turn  trajectory.  The  idea  of  gravity-turn  is 
that  the  vehicle  would  be  in  a  vertical  orientation,  and  then  nudged  over  a  small  amount 
so  that  the  center  of  mass  of  the  vehicle  is  no  longer  along  the  vertical  axis,  providing  a 
force  from  gravity  that  will  be  used  to  turn  the  vehicle. 

The  first  program,  fNITIALCOND,  is  an  m-file  with  a  section  noted  for  the  initial 
conditions.  This  file  would  provide  the  following: 

•  Initial  5  values  of  the  state  vector 

•  Time  step  and  final  time 

•  Prepare  the  output  file  to  receive  information 

•  Call  on  the  second  m-file  through  the  ODE45  command 

The  second  m-file  for  INITIALCOND  would  be  LAUNCHEOMS.  In 
EAUNCHEOMS,  the  equations  of  motion  are  written  in  terms  of  the  initial  state  vector. 
The  two  files  then  work  through  the  time  period,  which  in  this  case  is  designated  by  the 
bum  time  of  the  rocket  motor  of  82  seconds,  to  produce  a  final  state  vector  for  the  stage. 
The  information  is  output  to  a  text  file  at  each  time  step  so  that  the  data  can  later  be 
examined  and  plotted. 

At  the  end  of  INITIALCOND ’s  mn,  the  vehicle  has  expended  its  first  stage  and 
jettisons  the  dead  weight.  This  produces  a  discontinuity  in  the  change  of  mass  that 
dictates  the  need  to  write  separate  programs  for  each  of  the  stages.  The  second  program 
is  SECONDSTAGE  and  is  paired  with  SECONDEOMS.  The  program  SECONDSTAGE 
is  virtually  identical  to  INITIALCOND  except  that  it  calls  on  a  new  initial  mass,  and  then 
calls  on  the  last  line  of  INITIALCOND ’s  output  for  the  remainder  of  the  initial  state 
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vector.  Then  SECONDSTAGE  calls  SECONDEOMS,  which  again  runs  through  the 
equations  of  motion  provided  before.  SECONDEOMS  is  identical  to  LAUNCHEOMS 
except  in  the  variables  assigned  to  the  state  vector.  The  variables  were  named  differently 
to  provide  easier  access  for  data  manipulation.  At  the  end  of  SECONDSTAGE,  again 
there  is  a  discontinuity  in  mass  from  jettison  of  the  second  stage  dead  weight.  Also  note 
that  the  third  stage  has  a  different  thrust,  mass  flow  rate,  and  bum  time.  Therefore 
THIRDSTAGE  is  very  similar  to  its  predecessors  but  has  very  different  values.  It  also 
calls  on  the  last  vector  provided  by  SECONDSTAGE.  The  output  from  THIRDSTAGE 
is  then  provided  in  the  output  file,  where  a  five-column  matrix  is  given.  Each  stage  has 
it’s  own  text  file  for  easier  differentiation  of  the  data.  Of  primary  interest  in  the  results 
for  this  research  is  the  final  flight  path  angle.  Two  graphs  would  be  plotted  using  the  m- 
file,  PEOTFIEE.  The  first  graph  would  plot  altitude  vs.  downrange  distance,  which  in 
essence  is  a  modeling  of  the  actual  flight  path  of  the  vehicle  as  shown  in  Figure  4.  The 
second  graph  (Figure  5)  plots  altitude  vs.  flight  path  angle,  which  provides  a  graphical 
means  to  observe  the  behavior  of  the  vehicle’s  attitude.  With  an  initial  reference 
trajectory  provided,  we’re  ready  to  extrapolate  to  the  desired  point  of  a  zero  degree  flight 
path  angle  at  bum  out. 


20 


Altitude 

(meters) 


^  Altitude  vs  Distance  Stages  1  through  3 


(meters) 


Figure  5:  Altitude  vs  Downrange  Distance  for  Initial  yo  of  89,5° 
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Altitude 

(meters) 


^  Altitude  vs  Gamma  Stages  1  through  3 


Figure  6:  Altitude  vs  Gamma  y  for  yo  of  89,5" 


22 


3.3  Extrapolation  to  Zero  Flight  Path  Angle  at  Burn  Out 


The  goal  here  is  to  find  the  correct  initial  conditions  that  provide  for  a  zero  degree 
flight  path  angle  at  burnout.  To  do  this,  a  second  trajectory  is  found  using  the  same 
procedure  as  for  the  initial  trajectory  above.  However,  this  trajectory  begins  with  an 
angle  of  88.5”  as  it’s  initial  flight  path  angle  Yq.  After  running  the  m-files  through  with 
this  initial  condition,  a  second  trajectory  is  formed  as  shown  in  Figures  7  and  8  below. 


x10^ 


Distance  x10^ 

(meters) _ 


Figure  7:  Altitude  vs  Distance  for  yo  of  88,5" 
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Figure  8:  Altitude  vs  Flight  Path  Angle  jo  for  yo  of  88,5" 


As  shown  in  these  eharts,  the  final  y  is  not  equal  to  zero.  Here  another  m-file 
routine  ealled  NEWGAMMA  is  applied.  This  program  requires  values  to  be  input  into 
the  file  direetly.  The  initial  flight  path  angles,  yo,  and  the  final  angles,  yf,  from  both  runs 
are  input  here.  The  program  then  runs  a  simple  formula  to  find  the  next  yoi 


A/f 


rf. 


(3-1) 
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Where  7 is  the  last  initial  flight  path  angle,  Ayo  is  the  differenee  between  the  most 
reeent  two  initial  flight  path  angles,  Ayf  is  the  differenee  between  the  two  most  reeent 

final  flight  path  angles,  and  7 is  the  most  reeent  final  flight  path  angle  of  the  vehiele. 

This  last  term  would  be  a  Ay  of  the  differenee  between  the  given  and  desired  final  flight 
path  angle  exeept  our  final  desired  angle  is  zero.  Then  the  new  yo  is  input  into  the  next 
run  as  the  initial  flight  path  angle.  This  proeess  is  iterative  as  the  value  that  will  bring  yf 
to  zero  is  found.  The  entire  set  of  m-files  is  then  run  multiple  times  giving  the  results 
diseussed  in  the  following  ehapter. 
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IV.  RESULTS  AND  DISCUSSION 


4.1  Introduction 

This  chapter  presents  the  results  found  in  the  gravity-turn  trajectory  modeling  of 
the  Boeing  AirLaunch  based  vehiele.  The  only  payload  considered  in  this  researeh  was 
the  Spaee  Maneuver  Vehicle  (SMV)  that  is  currently  in  development  for  the  Air  Force. 
This  is  considered  to  be  approximately  3,000  kg  (7,500  lb)  of  payload  mass.  With  the 
vehicle  described  in  previous  ehapters  as  the  launch  vehicle,  the  gravity-turn  m-files  ean 
now  be  run. 


4.2  Initial  Results 

Since  the  engines  are  solid  rocket  motors,  they  are  not  capable  of  throttling  by 
design.  This  provides  that  there  is  no  capability  for  changing  thrust  during  the  burn 
sequences.  Intuition  should  prove  that  given  these  facts  this  method  has  a  fixed  altitude 
at  whieh  the  final  flight  path  angle  yf  goes  to  zero.  Using  the  proeedures  deseribed  in  the 
previous  ehapter,  the  first  runs  were  found  to  go  to  zero  at  the  state  deseribed  in  Table  3. 
This  run  is  referred  to  as  case  1  for  the  remainder  of  this  thesis. 
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Table  3:  Final  conditions  of  case  1 


Initial  Flight  Path  Angle  Yo 

87.51904" 

Altitude 

81433.5009  meters 

Downrange  Distance 

504090.589  meters 

Velocity 

7344.29656  m/s 

As  can  be  seen  in  Table  3,  the  final  altitude  is  only  81.4  kilometers.  It  is 
immediately  obvious  that  this  is  not  high  enough  to  be  in  low  earth  orbit.  For  this 
research,  an  orbit  of  250  kilometers  was  ehosen  as  the  target  altitude.  Several  ideas  were 
approaehed  to  alter  the  final  altitude.  Sinee  the  vehiele  is  powered  by  solid  roeket 
motors,  there  is  no  eapability  for  variation  of  thrust  or  bum  time.  One  possibility, 
therefore,  is  to  look  into  coasting  the  vehicle  between  stages  in  order  to  gain  a  higher 
altitude. 


4.3  Vehicle  Coasting 

In  order  to  reach  higher  orbit  altitudes,  coasting  the  vehicle  was  investigated.  This 
involves  waiting  between  stages  before  igniting  the  next  stage.  After  the  previous  stage 
has  burned  out,  it  will  be  jettisoned  and  the  vehicle  will  then  eontinue  on  for  a  specified 
amount  of  time  under  no  thmst.  This  is  referred  to  as  coasting  as  the  vehicle  is  under  no 
power.  By  coasting,  the  vehiele  will  eontinue  to  gain  altitude  while  conserving  fuel  and 
can  therefore  get  to  higher  altitudes.  To  insert  coasting,  m-files  were  ereated  that  ran  the 
equations  of  motion  without  thmst  or  mass  flow.  Similar  to  the  routines  deseribed  above. 
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the  initial  flight  path  angle  had  to  be  varied  until  a  final  flight  path  angle  of  zero  was 
attained.  Several  runs  were  attempted  while  varying  the  eoast  times  between  stage 
ignitions.  As  discussed  before,  drag  is  assumed  to  be  negligible,  which  could  cause 
differences  in  real  world  scenarios.  The  results  varied  as  presented  in  the  remainder  of 
this  section. 

The  first  cases  attempted  were  made  while  keeping  the  coast  times  consistent  with 
each  other.  In  case  2,  only  a  five  second  coast  was  considered.  As  seen  in  Figure  9  and 
Table  4  this  showed  that  only  a  small  change  was  attained.  Obviously  a  much  larger 
coast  time  is  needed. 


Table  4:  Final  conditions  of  case  2 


Initial  Flight  Path  Angle  Yo 

87.88“ 

Altitude 

86441.6828  meters 

Downrange  Distance 

505497.292  meters 

Velocity 

7263.80507  m/s 
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Figure  9:  Launch  trajectory  for  case  2 


From  case  2,  a  jump  to  a  coast  of  thirty-five  seconds  between  stages  was 
performed  in  case  3  to  see  how  large  a  difference  could  be  found.  As  seen  in  Figure  10 
and  Table  5,  a  much  higher  altitude  is  reached  than  the  previous  runs.  However,  an 
altitude  of  121  km  is  still  fairly  low  for  an  orbit. 


29 


Figure  10:  Launch  trajectory  for  case  3 


Table  5:  Final  conditions  for  case  3 


Initial  Flight  Path  Angle  Yo 

89.57“ 

Altitude 

121398.648  meters 

Downrange  Distanee 

464387.699  meters 

Veloeity 

6644.82810  m/s 
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In  case  4,  a  coast  time  of  sixty  seconds  between  stages  was  investigated.  This 
however  had  some  serious  problems  as  seen  in  Figure  1 1  and  Table  6.  It  was  found  that 
the  flight  path  did  not  converge  to  zero  for  this  coast  time,  and  the  vehicle  would  impact 
the  surface  during  the  second  stage  burn.  This  is  obviously  a  bad  thing  to  happen  to  the 
launch  vehicle. 
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Figure  11:  Launch  trajectory  for  case  4 
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Table  6:  Final  conditions  after  case  4 


Initial  Flight  Path  Angle  Yo 

89.9° 

Altitude 

-70583.0107  meters 

Downrange  Distance 

26483.1574  meters 

Velocity 

3.24739556  m/s 

After  several  experimental  runs,  it  was  determined  that  a  eoast  time  larger  than 
approximately  50  seeonds  during  the  first  eoast  eaused  the  vehiele’s  flight  path  angle  to 
degrade  too  quiekly.  The  deeision  was  then  made  to  confine  the  first  coast  time  to  50 
seconds  or  less.  Then  the  second  coast  times  were  varied  to  attain  higher  altitudes.  In 
case  5,  the  first  coast  stage  was  set  to  40  seconds  and  then  run  with  second  coast  times  of 
60  seconds,  90  seconds  and  120  seconds  which  are  represented  in  Figures  12,  13,  and  14, 
respectively.  These  figures  give  a  representation  of  possible  altitudes  to  reach,  as  all 
altitudes  between  them  are  capable  of  being  achieved  through  modification  of  the  coast 
times.  Also  note  the  large  difference  in  altitudes  between  the  different  coast  times.  This 
gives  evidence  that  the  vehicle  has  feasibility  to  reach  more  reasonable  orbital  altitudes. 
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Figure  12:  Launch  trajectory  of  case  5  with  second  stage  coast  of  60  seconds 
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Figure  13:  Launch  trajectory  of  case  5  with  a  second  stage  coast  of  90  seconds 
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Figure  14:  Launch  trajectory  of  case  5  with  a  second  stage  coast  of  120  seconds 


The  highest  altitude  achieved  through  the  case  5  trajectories  was  approximately 
237  km.  This  is  much  closer  to  the  more  desirable  orbits.  Attempts  to  reach  higher  orbits 
with  case  5’s  first  coast  were  unproductive.  So  to  achieve  slightly  higher  altitudes  in  case 
6,  the  first  coast  time  was  extended  to  50  seconds  and  the  second  coast  stage  time  was  set 
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at  120  seconds.  This  second  stage  coast  time  gave  the  highest  altitude  previously,  and  is 
therefore  a  logical  starting  place.  The  results  of  case  6  are  shown  in  Figure  15  and  Table 


7  below. 


5 

X10 


Figure  15:  Launch  trajectory  of  case  6 


Table  7:  Final  conditions  of  case  6 


Initial  Flight  Path  Angle  Yo 

89.98” 

Altitude 

228092.170  meters 

Velocity 

5450.60642  m/s 
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As  seen  in  Figure  14  and  Table  7,  the  altitude  reaehed  is  actually  less  than  that  of 
the  highest  achievable  in  case  5.  After  several  iterations,  a  second  coast  stage  of  160 
seconds  chosen  for  case  7.  The  trajectory  and  final  conditions  of  case  7  are  given  in 
Figure  16  and  Table  8,  respectively. 
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Figure  16:  Launch  trajectory  with  case  7 
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Table  8:  Final  conditions  after  case  7 


Initial  Flight  Path  Angle  Yo 

89.9975“ 

Altitude 

262427.963  meters 

Velocity 

4584.99328  m/s 

As  seen  in  Table  8,  the  vehicle  is  finally  above  250  km  in  altitude.  Now  that 
decent  altitudes  are  capable  of  being  reached,  the  vehicle  was  then  set  up  to  achieve  a 
certain  altitude. 

For  case  8,  a  goal  was  set  for  a  250  km  orbit.  A  trajectory  was  then  found  that 
would  approach  250  km  as  the  flight  path  angle  came  to  zero  and  then  the  vehicle  was 
allowed  to  keep  a  zero  flight  path  angle  and  accelerate  in  the  orbit.  The  results  of  case  8 
produce  the  trajectory  and  final  conditions  shown  in  Figure  17  and  Table  9  below. 
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Downrange  Distance  (meters) 


Figure  17:  Launch  trajectory  of  case  8 


Table  9:  Final  conditions  of  case  8 


Of  interest  to  the  researeh  at  this  point  is  the  velocity  needed  to  maintain  a 
circular  orbit.  This  velocity  is  governed  by  the  equation: 


where  p  is  the  gravitational  parameter  of  the  local  earth  system,  and  is  given  as 
3.98601x10^  kmVs^,  and  r  =  (Rg  +  H),  with  Rg  given  as  6378.135  km.  Assuming  an 
altitude  (H)  of  250  km,  it  is  found  that  the  velocity  required  to  maintain  this  circular  orbit 
is  approximately  7.754  km/s.  Given  the  final  velocity  shown  in  Table  9  of  5.945  km/s,  it 
is  obvious  that  the  vehicle  will  not  remain  in  this  orbit.  So  an  addition  of  velocity  is 
required  to  maintain  the  desired  orbit.  As  previously  stated,  the  payload  of  the  Boeing 
AirLaunch  System  for  this  research  is  assumed  to  be  the  Space  Maneuvering  Vehicle 
(SMV).  The  SMV  is  still  in  development,  however  by  design  it  will  be  capable  of  3200 
m/s  of  additional  velocity  commonly  referred  to  as  AV  (“Space  Maneuver  Vehicle  Fact 
Sheet”:  2004).  For  the  purpose  of  reaching  the  desired  orbital  velocity,  it  is  assumed  that 
the  SMV  will  be  used  to  produce  the  required  AV  of  1809  m/s,  which  is  well  within  its 
capabilities.  Assuming  the  SMV  accelerates  at  the  same  rate  that  the  third  stage  did, 
about  100  m/s  ,  this  will  add  approximately  18  seconds  to  the  launch  time.  This  gives  the 
time  from  ignition  to  on  orbit  as  approximately  429  seconds  or  7.15  minutes. 
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4.4  Hohmann  Transfer  Results 


Another  possibility  for  orbital  insertion  is  a  transfer  orbit  into  the  desired  plane. 
Sinee  this  vehiele  has  a  limited  amount  of  fuel  available  to  it,  a  Hohmann  transfer  is  the 
best  ehoice  to  make  the  transfer  when  time  to  orbit  is  not  an  issue.  As  discussed  in 
Chapter  2,  the  Hohmann  transfer  consists  of  two  changes  in  velocity  to  first  put  the 
vehicle  onto  an  elliptical  orbit  and  then  another  to  place  the  orbiter  into  the  desired  orbit. 

For  the  first  consideration,  case  1  was  taken  as  the  launch  trajectory  and  a  250  km 
orbit  again  was  chosen  as  the  target.  This  investigation  is  here  dubbed  case  9.  First,  the 
vehicle  at  the  end  of  case  1  is  at  approximately  81.43  km  and  traveling  at  7.344  km/s.  In 
order  to  take  advantage  of  equations  2-14  and  2-15,  the  vehicle  must  first  be  accelerated 
to  the  correct  velocity  to  maintain  the  given  orbit.  This  can  be  found  using  equation  4-1 
and  determining  the  difference  in  velocities.  Equations  2-14  and  2-15  can  then  be  used  to 
determine  the  AV’s  required  for  the  transfer,  and  equation  2-16  is  used  to  find  the  time  of 
flight  of  the  vehicle.  The  results  are  shown  below  in  Table  10. 


Table  10:  Case  9  results  for  Hohmann  transfer 


AV  required  for  inner  orbit 

0.51 1  km/s 

AV 1  for  the  Hohmann  transfer 

0.0504  km/s 

AV2  for  the  Hohmann  transfer 

0.0501  km/s 

Total  AV  required 

0.61 15  km/s 

Time  of  flight.  At 

2634.09  s  or  43.9  minutes 
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With  a  total  AV  required  of  0.61 15  km/s,  this  maneuver  is  well  within  the  range 
of  the  SMV’s  eapability  of  a  3.2  km/s  AV.  This  shows  an  exeellent  possibility  for 
launehing  to  orbit. 

In  ease  10,  the  vehiele  is  assumed  to  be  at  the  end  of  ease  6.  At  an  altitude  of 
approximately  228.092  km  and  veloeity  of  5.45  km/s,  the  AV’s  required  to  achieve  a  250 
km  orbit  are  shown  in  Table  1 1 . 


Table  11;  Case  10  results  for  Hohmann  transfer 


AV  required  for  inner  orbit 

2.3177  km/s 

AV  1  for  the  Hohmann  transfer 

0.00643  km/s 

AV2  for  the  Hohmann  transfer 

0.00642  km/s 

Total  AV  required 

2.3305  km/s 

Time  of  flight.  At 

2678.49  s  or  44.6  min 

The  results  of  case  10  are  still  within  the  range  of  the  available  AV  of  the  SMV, 
however,  case  9  is  a  much  more  efficient  profile.  Interestingly  this  poses  an  unexpected 
scenario.  The  initial  launch  trajectory  posed  in  case  1  barely  reaches  the  altitude  of  80 
km  but  is  far  more  efficient  using  a  Hohmann  transfer  than  coasting  is,  as  displayed  in 
case  8.  Also  of  interest  here  is  that  the  higher  trajectories  reached  by  coasting  have  a 
higher  total  AV  than  the  trajectory  that  has  no  coasting.  This  is  primarily  due  to  the  fact 
that  as  the  vehicle  coasts  it  is  slowing  down  due  to  gravity  losses  while  the  launch 
trajectory  without  coasting  provides  a  velocity  very  near  the  orbital  velocity  needed. 
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The  results  shown  in  this  ehapter  demonstrate  that  the  vehiele  represented  in  this 
research  is  capable  of  reaching  a  prescribed  altitude  through  a  gravity-turn  trajectory  with 
supplementation.  Surprising  to  this  researcher,  the  initial  launch  augmented  with  a 
Hohmann  transfer  is  the  most  fuel  efficient  of  the  cases  investigated.  Using  the 
capability  of  an  air  launch  platform,  this  vehicle  is  capable  of  intercepting  the  proper 
orbit  of  a  target  and  launching  at  such  a  time  that  at  engines  cut  off,  the  target  is 
alongside.  For  a  fuel-efficient  attempt,  the  vehicle  could  rendezvous  at  the  end  of  the 
Hohmann  transfer,  and  for  a  faster  intercept  the  vehicle  can  use  coasting  to  achieve  the 
proper  orbital  altitude  and  use  the  SMV  to  accelerate  into  the  orbital  plane,  as  was  shown 
in  case  8.  The  following  chapter  discusses  this  as  well  as  recommendations  for  future 
investigation. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


5.1  Conclusions 

The  goal  of  this  research  was  to  evaluate  the  feasibility  of  a  coplanar,  gravity-turn 
launch  trajectory  using  a  vehicle  based  heavily  upon  the  Boeing  AirLaunch  System.  The 
results  from  the  previous  chapter  demonstrate  that  the  vehicle  described  in  this  research  is 
capable  of  reaching  reasonable  altitudes  in  low  earth  orbit  with  modifications  to  its 
trajectory  through  coasting  and  transfer  orbits.  Since  the  payload  of  the  vehicle  can  be 
considered  in  essence  a  fourth  stage,  there  are  a  considerable  number  of  options  that  can 
be  used  to  achieve  a  desired  orbit.  This  research  limited  these  options  to  coasting  and 
Hohmann  transfers. 

Since  the  vehicle  can  be  positioned  literally  anywhere  in  the  world,  the  launch 
scenarios  investigated  are  quite  possible.  With  actual  information  on  the  third  stage  and 
the  Space  Maneuver  Vehicle  (SMV),  calculations  of  a  direct  to  rendezvous  launch  is 
possible.  This  would  simply  require  determining  the  correct  launch  time  as  dictated  by 
the  target  vehicle.  The  mission  could  be  completed  in  however  long  it  would  take  the 
Boeing  747  launch  vehicle  to  arrive  at  the  proper  location. 

This  research  provides  an  interesting  and  insightful  look  at  gravity-turn  launch 
trajectories  and  the  use  of  solid  rocket  motors  in  a  responsive  launch  situation.  It  shows 
that  with  coasting  and  transfer  orbits,  most  low  earth  orbits  are  attainable,  and  direct  to 
rendezvous  is  feasible.  As  an  interesting  development,  it  also  demonstrated  that  the  most 
efficient  launch  was  to  a  lower  altitude  followed  by  a  Hohmann  transfer  to  the  desired 
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altitude.  The  tradeoff  here  of  eourse  is  time  to  rendezvous  at  the  eorreet  orbit.  Case  8 
achieves  engines  off  on  orbit  at  just  over  7  minutes,  where  as  case  9  reaches  the  target 
orbit  at  approximately  48  minutes  after  engine  ignition. 


5.2  Recommendations 

The  most  obvious  next  step  for  this  system  is  to  find  actual  specifications  for  the 
third  stage  and  SMV.  This  may  greatly  affect  the  system  if  actual  specifications  are 
much  greater  or  less  than  the  assumed  values.  Also  of  interest  is  to  add  drag  into  the 
equation.  The  introduction  of  drag  will  primarily  affect  the  first  stage  and  possibly  the 
first  coasting  stage.  Another  point  of  interest  would  be  to  reproduce  the  equations  of 
motion  in  the  three  dimensional  earth  centered  reference  frame.  With  the  three- 
dimensional  models,  a  launch  to  rendezvous  can  be  examined  in  depth  for  the  purposes  of 
mission  planning. 

As  the  field  of  space  launch  continues  to  evolve,  the  desirable  capabilities  of 
launch  on  demand  and  direct  to  rendezvous  launch  will  become  more  and  more  feasible 
in  reality.  Platforms  like  the  Boeing  AirLaunch  System  are  not  too  far  into  the  future. 
Orbital  already  operates  a  much  smaller  payload  class  on  a  similar  platform  called 
Pegasus.  The  results  of  this  thesis  offer  an  insight  into  the  possible  next  evolution  of  this 
developing  technology.  The  ability  of  launching  directly  to  rendezvous  has  significant 
advantages  in  civilian  and  military  aspects  as  well  as  peacetime  and  wartime  missions. 
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The  portability  of  an  air  launch  system  allows  the  concept  of  reaching  virtually  any 
orbital  plane  desirable.  The  research  presented  here  poses  some  interesting  possibilities 
for  this  type  of  technology  in  the  future. 
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APPENDIX  A:  Code  Summary 


FIRSTCOASTEOM 

TYPE:  Subroutine  m-file 

PURPOSE:  This  routine  provides  the  equations  of  motion  to  iterate  through  the  first 

coasting  stage  conditions.  This  file  provides  the  changes  in  the  state 
variables  per  time. 

INPUTS:  State  vector 

OUTPUTS:  State  vector 

CALLS:  None 

AUTHOR:  David  W.  Callaway 

FIRSTSTAGECOAST 

TYPE:  Main  Program  m-file 

PURPOSE:  This  program  solves  a  boundary  value  problem  for  the  gravity-turn 

trajectory  of  the  first  coasting  stage. 

INPUTS:  Final  state  vector  from  first  stage  bum  and  new  initial  mass 

OUTPUTS:  Final  state  vector  after  coasting 
CALLS:  FIRSTCOASTEOM 

INITIALCOND 

AUTHOR:  David  W.  Callaway 
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INITIALCOND 


TYPE:  Main  Program  m-file 

PURPOSE:  This  program  solves  a  boundary  value  problem  for  the  gravity-turn 

trajeetory  of  the  first  stage  burning. 

INPUTS:  Initial  state  veetor 

OUTPUTS:  Pinal  state  vector  for  first  stage 

CAPES:  EAUNCHEOMS 

AUTHOR:  David  W.  Callaway 

EAUNCHEOMS 

TYPE:  Subroutine  m-file 

PURPOSE:  This  routine  provides  the  equations  of  motion  to  iterate  through  the  first 

stage  bum  conditions.  This  file  provides  the  changes  in  the  state  variables 
per  time. 

INPUTS:  State  vector 

OUTPUTS:  State  vector 

CALES:  None 

AUTHOR:  David  W.  Callaway 
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NEWGAMMA 


TYPE: 

PURPOSE: 

INPUTS: 

OUTPUTS: 

GALES: 

AUTHOR: 

PLOTEILE 

TYPE: 

PURPOSE: 

INPUTS: 

OUTPUTS: 

CALLS: 


AUTHOR: 


Main  Program  m-file 

This  program  determines  the  ehange  in  final  flight  path  angle  per  ehange 
in  initial  flight  path  angle.  It  then  finds  the  new  initial  flight  path  angle  for 
the  next  iteration. 

Initial  flight  path  angles  from  the  two  previous  iterations 
Pinal  flight  path  angles  from  the  two  previous  iterations 

Initial  flight  path  angle 

None 

David  W.  Callaway 


Main  Program  m-file 

This  program  ealls  the  iterated  state  veetors  from  all  stages  of  the  launeh 
and  plots  them  for  data  analysis. 

All  Iterations  of  state  veetors 

Downrange  distanee  vs.  Altitude  plot 
Plight  path  angle  vs.  Altitude  plot 

output  from: 

INITIALCOND 

PIRSTSTAGECOAST 

SECONDSTAGE 

SECONDSTAGECOAST 

THIRDSTAGE 

David  W.  Callaway 
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SECONDCOASTEOMS 


TYPE: 

PURPOSE: 

INPUTS: 

OUTPUTS: 

GALES: 

AUTHOR: 


Subroutine  m-file 

This  routine  provides  the  equations  of  motion  to  iterate  through  the 
second  coasting  stage  conditions.  This  file  provides  the  changes  in  the 
state  variables  per  time. 

State  vector 

State  vector 

None 

David  W.  Callaway 


SECONDEOMS 


TYPE: 

PURPOSE: 

INPUTS: 

OUTPUTS: 

CALLS: 

AUTHOR: 


Subroutine  m-file 

This  routine  provides  the  equations  of  motion  to  iterate  through  the 
second  stage  burn  conditions.  This  file  provides  the  changes  in  the  state 
variables  per  time. 

State  vector 

State  vector 

None 

David  W.  Callaway 
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SECONDSTAGE 


TYPE: 

Main  Program  m-file 

PURPOSE: 

This  program  solves  a  boundary  value  problem  for  the  gravity-turn 
trajectory  of  the  second  stage  burning. 

INPUTS: 

Pinal  state  vector  of  either  first  stage  or  first  coasting  stage,  dependent 
on  case. 

OUTPUTS: 

Pinal  state  vector  for  second  stage 

CAEES: 

SECONDEOMS 

AUTHOR: 

David  W.  Callaway 

SECONDSTAGECOAST 


TYPE: 

PURPOSE: 

INPUTS: 

OUTPUTS: 

CAEES: 

AUTHOR: 


Main  Program  m-file 

This  program  solves  a  boundary  value  problem  for  the  gravity-turn 
trajeetory  of  the  second  coasting  stage. 

Pinal  state  vector  of  second  stage 

Pinal  state  vector  for  second  coasting  stage 

SECONDCOASTEOMS 

David  W.  Callaway 
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THIRDEOMS 


TYPE: 

Subroutine  m-file 

PURPOSE: 

This  routine  provides  the  equations  of  motion  to  iterate  through  the  third 
stage  bum  conditions.  This  file  provides  the  changes  in  the  state  variables 
per  time. 

INPUTS: 

State  vector 

OUTPUTS: 

State  vector 

GALES: 

None 

AUTHOR: 

David  W.  Callaway 

THIRDSTAGE 

TYPE: 

Main  Program  m-file 

PURPOSE: 

This  program  solves  a  boundary  value  problem  for  the  gravity-turn 
trajectory  of  the  third  stage  burning. 

INPUTS: 

Pinal  state  vector  of  second  stage  or  second  coasting  stage,  dependent  on 
case 

OUTPUTS: 

Pinal  state  vector  for  third  stage 

CALLS: 

THIRDEOMS 

AUTHOR: 

David  W.  Callaway 
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APPENDIX  B:  M-Files 


The  m-files  presented  here  are  from  ease  7.  It  should  be  noted  that  they  are 
provided  here  in  alphabetieal  order;  however,  they  must  be  run  from  the  same  file  in  the 
order  INITIALCOND,  FIRSTSTAGECOAST,  SECONDSTAGE, 
SECONDSTAGECOAST,  THIRDSTAGE,  and  PEOTEIEE. 


EIRSTCOASTEQM 


%  eoms  for  firststage  eoast 
%ET  Callaway,  Thesis 
%Trajeetory  model 

funetion  [yedot]=firsteoasteom(t,ye) 
global  mdot  T  Re  ge  tstep 

%EOM's  for  2-D  launeh 

%yO=[mnot  Knot  Xnot  Vnot  gammanot}; 

%ohange  in  mass  mdot 
yodot(l,l)=0; 

%ohange  in  altitude  Hdot 
yedot(2,l)=ye(4)*sin(ye(5)); 

%change  in  horizontal  distance  Xdot 
ycdot(3 , 1  )=yc(4)*cos(yc(5)); 

%change  in  Velocity 

ycdot(4,l)=(l/yc(l))*(T-(yc(l)*ge-yc(l)*((yc(4)*cos(yc(5)))^2)/(Re+yc(2)))*sin(yc(5))); 
%change  in  gamma 

ycdot(5,l)=-(l/yc(l))*(l/yc(4))*(yc(l)*ge- 

yc(l)*((yc(4)*cos(yc(5)))^2)/(Re+yc(2)))*cos(yc(5)); 
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FIRSTSTAGECQAST 
%  coast  stage  1 
%LT  Callaway,  Thesis 
%Trajeetory  model 
%Intial  Conditions: 

elose  all 

fde=fopen('C:\thesis\withooast\firststageooast.txtVw+'); 
fprintf(lile,'mass  \t  \t  \t  altitude  \t  \t  \t  X  \t  \t  \t  V  \t  \t  \t  gamma  \n'); 

global  mdot  T  Re  ge  tstep 

%Times 

tstep=l;  %see 

tfmal=50;  %sec 


%  Constants - 

%Thrust  in  N 

%mass  flow  is  eonsidered  eonstant  (kg/s); 
%Radius  of  earth  in  meters; 

%gravity  g; 


T=0; 

mdot=0; 

Re=6378135; 

ge=9.81; 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


%  Initial  Conditions 


mnot=83388. 874893;  %total  mass  -  first  stage 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


%  yyyO=[mnot  Knot  Xnot  Vnot  gammanot]; 
yeO=[mnot  y(83,2)  y(83,3)  y(83,4)  y(83,5)]; 

%  Integrate  the  equations  of  motion 

options  =  odeset('RelTor,le-7,AbsTor,le-10*ones(l,5)); 

[t,yo]  =  ODE45('firstooasteom',0:tstep:tfmal,yoO, options); 

out=[yo(:,2)  yo(:,3)  yc(:,4)  yc(:,5)]; 
fprmtf(fde,'%2.8e  \t  %2.8e  \t  %2.8e  \t  %2.8e  \n',out'); 

%  fprmtf(fde,’%2.8e  \t  %2.8e  \t  %2.8e  \t  %2.8e  \t  %2.8e  \n’,yc’); 
folose(fde); 
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INITIALCOND 
%LT  Callaway,  Thesis 
%Trajectory  model 
%Intial  Conditions: 

elear 
elose  all 

file=fopen('C:\thesis\withooast\firststage.txtVw+'); 
fprintf(lile,'mass  \t  \t  \t  altitude  \t  \t  \t  X  \t  \t  \t  V  \t  \t  \t  gamma  \n'); 

global  mdot  T  Re  ge  tstep 

%Times 
tstep=l; 
tfinal=82; 

%  Constants - 

1=1650245.73705; 
mdot=592. 64056; 

Re=6378135; 
ge=9.81; 


%Thrust  in  N 

%mass  flow  is  eonsidered  eonstant  (kg/s); 
%Radius  of  earth  in  meters; 

%gravity  g; 


%see 

%sec 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


%  Initial  Conditions 


mnot=136077.711; 

gammanot=89.9975*pi/180; 

Vnot=134.112; 

Hnot=6096; 

Xnot=0; 


%total  mass  in  kg 
%gamma  must  be  in  radians; 

%Initial  veloeity  300  mph,  but  in  meters/see: 
%Initial  Altitude...  20kft...m  meters 
%Initial  X  position 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


y0=[mnot  Knot  Xnot  Vnot  gammanot]; 


%  Integrate  the  equations  of  motion 

options  =  odeset('RelTor,le-7,'AbsTor,le-10*ones(l,5)); 

[t,y]  =  ODE45('launoheoms',0:tstep:tfmal,yO, options); 


out=[y(:,2)  y(:,3)  y(:,4)  y(:,5)]; 

fprmtf(fde,'%2.8e  \t  %2.8e  \t  %2.8e  \t  %2.8e  \n',out'); 

%fprintf(fde,’%2.8e  \t  %2.8e  \t  %2.8e  \t  %2.8e  \t  %2.8e  \n',y’); 

folose(fde); 
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LAUNCHEOMS 
%  Launcheoms  for  first  stage 
%  LT  Callaway,  Thesis 
%Trajeetory  model 

function  [ydot]=launcheoms(t,y) 
global  mdot  T  Re  ge  tstep 

%EOM's  for  2-D  launch 

%yO=[mnot  Knot  Xnot  Vnot  gammanot}; 

%change  in  mass  mdot 
ydot(l,l)=-592.64056; 

%change  in  altitude  Hdot 
ydot(2,l)=y(4)*sm(y(5)); 

%change  in  horizontal  distance  Xdot 
y  dot(3 , 1  )=y(4)  *  cos(y(5)); 

%change  in  Velocity 

ydot(4,l)=(l/y(l))*(T-(y(l)*ge-y(l)*((y(4)*cos(y(5)))^2)/(Re+y(2)))*sm(y(5))); 
%change  in  gamma 

ydot(5,l)=-(l/y(l))*(l/y(4))*(y(l)*ge-y(l)*((y(4)*cos(y(5)))^2)/(Re+y(2)))*cos(y(5)); 


NEWGAMMA 

%dgamma/dgammanot  program 
%dgamma=  old  gamma  -  new  gamma 
gammao=5.25671725e-001; 
gamman=1.15971452e-002; 
dgamma=gammao-gamman  %radians 

%  dgammanot=  |new  gamma  -  previous  gamma| 

gammanotnew=89 . 7 8 ; 

gammanotold=89.9; 

dgammanot=gammanotold-gammanotnew  %degrees 
deltagammanot=(dgammanot)  *  (gamman)/ (dgamma)  %degrees 
newgammanot=gammanotnew-deltagammanot 
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PLOTFILE 
hold  on; 

%plotting  altitude  vs  X 

plot(y(:,3),y(:,2)); 

plot(yc(:,3),yc(:,2)); 

plot(yy(:,3),yy(:,2)); 

plot(ycc(:,3),ycc(:,2)); 

plot(yyy(:,3),yyy(:,2)); 

%plotting  altitude  vs  gamma 
%  plot(y(:,5),y(:,2)); 

%  plot(ye(:,5),ye(:,2)); 

%  plot(yy(:,5),yy(:,2)); 

%  plot(yee(:,5),yec(:,2)); 

%  plot(yyy(:,5),yyy(:,2)); 


SECONDCOASTEOMS 
%  eoms  for  seeond  stage  coast 

function  [yccdot]=secondcoasteoms(t,ycc) 
global  mdot  T  Re  ge  tstep 

%EOM's  for  2-D  launch 

%yO=[mnot  Knot  Xnot  Vnot  gammanot}; 

%change  in  mass  mdot 
yccdot(l,l)=0; 

%change  in  altitude  Hdot 
yccdot(2,l)=ycc(4)*sin(ycc(5)); 

%change  in  horizontal  distance  Xdot 
yccdot(3,l)=ycc(4)*cos(ycc(5)); 

%change  in  Velocity 

yccdot(4,l)=(l/ycc(l))*(T-(ycc(l)*ge- 

ycc(l)*((ycc(4)*cos(ycc(5)))^2)/(Re+ycc(2)))*sin(ycc(5))); 

%change  in  gamma 

yccdot(5,l)=-(l/ycc(l))*(l/ycc(4))*(ycc(l)*ge- 

ycc(l)*((ycc(4)*cos(ycc(5)))^2)/(Re+ycc(2)))*cos(ycc(5)); 
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SECONDEOMS 
%  Eauncheoms  for  second  stage 

function  [yydot]=secondeoms(t,yy) 
global  mdot  T  Re  ge  tstep 

%EOM's  for  2-D  launch 

%yO=[mnot  Knot  Xnot  Vnot  gammanot}; 

%change  in  mass  mdot 
yydot(l,l)=-592.64056; 

%change  in  altitude  Hdot 
yydot(2,  l)=yy(4)*sin(yy(5)); 

%change  in  horizontal  distance  Xdot 
yydot(3 , 1  )=yy  (4)  *  cos(yy(5)); 

%change  in  Velocity 

yydot(4, 1  )=(  1  /yy  ( 1 ))  *  (T-(yy  ( l)*ge- 

yy(l)*((yy(4)*cos(yy(5)))^2)/(Re+yy(2)))*sin(yy(5))); 

%change  in  gamma 

yydot(5,l)=-(l/yy(l))*(l/yy(4))*(yy(l)*ge- 

yy(l)*((yy(4)*cos(yy(5)))^2)/(Re+yy(2)))*cos(yy(5)); 
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SECONDSTAGE 

%  SECOND  STAGE  INIITIAL  CONDITIONS 
%Intial  Conditions: 

close  all 

file=fopen('C:\thesis\withcoast\secondstage.txtVw+'); 
fprintf(fde,'mass  \t  \t  \t  altitude  \t  \t  \t  X  \t  \t  \t  V  \t  \t  \t  gamma  \n'); 

global  mdot  T  Re  ge  tstep 

%Times 
tstep=l; 
tlinal=82; 

%  Constants - 

T=1650245.73705; 
mdot=592. 64056; 

Re=6378135; 
ge=9.81; 


%Thrust  in  N 

%mass  flow  is  considered  constant  (kg/s); 
%Radius  of  earth  in  meters; 

%gravity  g; 


%sec 

%sec 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


%  Initial  Conditions 


mnot=83388. 874893; 

%  gammanot=l. 36707253; 
%  Vnot=564.206512; 

%  Hnot=30732.3097; 
%Xnot=2685.21896; 


%total  mass  -  first  stage  mass  in  kg 
%gamma  must  be  in  radians; 

%Velocity  from  end  of  1st  stage,  but  in  meters/sec: 
%Altitude  from  previous  run 
%Initial  X  position 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


%  yy0=[mnot  Knot  Xnot  Vnot  gammanot]; 
yy0=[mnot  yc(5 1 ,2)  yc(5 1 ,3)  yc(5 1 ,4)  yc(5 1,5)]; 

%  Integrate  the  equations  of  motion 

options  =  odeset('RelTor,le-7,AbsTor,le-10*ones(l,5)); 

[t,yy]  =  ODE45('secondeoms',0:tstep:tfmal,yyO, options); 

out=[yy(:,2)  yy(:,3)  yy(:,4)  yy(:,5)]; 
fprintf(fde,'%2.8e  \t  %2.8e  \t  %2.8e  \t  %2.8e  \n',out'); 
%fprintf(fde,'%2.8e  \t  %2.8e  \t  %2.8e  \t  %2.8e  \t  %2.8e  \n',yy'); 
fclose(fde); 
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SECONDSTAGECOAST 


%  coast  stage  2 
%LT  Callaway,  Thesis 
%Trajectory  model 
%Intial  Conditions: 

elose  all 

file=fopen('C:\thesis\withooast\seoondstageooast.txtVw+'); 
fprintf(lile,'mass  \t  \t  \t  altitude  \t  \t  \t  X  \t  \t  \t  V  \t  \t  \t  gamma  \n'); 

global  mdot  T  Re  ge  tstep 

%Times 

tstep=l;  %see 

tfmal=160;  %see 


%  Constants - 

%Thrust  in  N 

%mass  flow  is  considered  eonstant  (kg/s); 
%Radius  of  earth  in  meters; 

%gravity  g; 


T=0; 

mdot=0; 

Re=6378135; 

ge=9.81; 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


%  Initial  Conditions 


mnot=30700. 027786;  %total  mass  -  stages  1  and  2 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


%  yyy0=[mnot  Knot  Xnot  Vnot  gammanot]; 
yee0=[mnot  yy(83,2)  yy(83,3)  yy(83,4)  yy(83,5)]; 


%  Integrate  the  equations  of  motion 

options  =  odeset('RelTor,le-7,AbsTor,le-10*ones(l,5)); 

[t,yoc]  =  ODE45('secondcoasteoms',0:tstep:tfmal,ycoO, options); 


out=[yco(:,2)  yoc(:,3)  ycc(:,4)  ycc(:,5)]; 
fprmtf(fde,'%2.8e  \t  %2.8e  \t  %2.8e  \t  %2.8e  \n',ouf); 

%  fprmtf(fde,'%2.8e  \t  %2.8e  \t  %2.8e  \t  %2.8e  \t  %2.8e  \n',ycc'); 
fclose(fde); 
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THIRDEOMS 


%  Launcheoms  for  third  stage 
%  Corrected  for  mdot=50%  of  mdot  stages  1  and  2 

function  [yyydot]=thirdeoms(t,yyy) 
global  mdot  T  Re  ge  tstep 

%EOM's  for  2-D  launch 

%yO=[mnot  Knot  Xnot  Vnot  gammanot}; 

%change  in  mass  mdot 
yyydot(l,l)=-296.32028; 

%change  in  altitude  Hdot 
yyydot(2,  l)=yyy(4)*sin(yyy(5)); 

%change  in  horizontal  distance  Xdot 
yyydot(3,l)=yyy(4)*cos(yyy(5)); 

%change  in  Velocity 

yyydot(4,l)=(l/yyy(l))*(T-(yyy(l)*ge- 

yyy(l)*((yyy(4)*cos(yyy(5)))^2)/(Re+yyy(2)))*sin(yyy(5))); 

%change  in  gamma 

yyydot(5,l)=-(l/yyy(l))*(l/yyy(4))*(yyy(l)*ge- 

yyy(l)*((yyy(4)*cos(yyy(5)))^2)/(Re+yyy(2)))*cos(yyy(5)); 
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THIRDSTAGE 

%  Third  STAGE  INIITIAL  CONDITIONS 
%Intial  Conditions: 

close  all 

file=fopen('C:\thesis\withcoast\thirdstage.txtVw+'); 
fprintf(fde,'mass  \t  \t  \t  altitude  \t  \t  \t  X  \t  \t  \t  V  \t  \t  \t  gamma  \n'); 

global  mdot  T  Re  ge  tstep 

%Times 

tstep=l;  %sec 

tfmal=84;  %see  rounded  down  from  84.9  seconds 


%  Constants - 

%Thrust  in  N 

%mass  flow  is  considered  constant  (kg/s); 
%Radius  of  earth  in  meters; 

%gravity  g; 


T=813932.5451; 
mdot=592. 64056; 
Re=6378135; 
ge=9.81; 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


%  Initial  Conditions 


mnot=30700.027786; 

%  gammanot=  1. 13822 111; 
%Vnot=224 1.82678; 

%  Hnot=124507.332; 
%Xnot=375 16.3229; 


%total  mass  -  first  stage  mass  in  kg 
%gamma  must  be  in  radians; 

%Velocity  from  end  of  2nd  stage,  in  meters/sec: 
%Altitude  from  previous  run 
%Initial  X  position  third  stage 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


%  yyy0=[mnot  Knot  Xnot  Vnot  gammanot]; 
yyy0=[mnot  ycc(161,2)  yec(161,3)  yec(161,4)  ycc(161,5)]; 

%  Integrate  the  equations  of  motion 

options  =  odeset('RelTor,le-7,AbsTor,le-10*ones(l,5)); 

[t,yyy]  =  ODE45('thirdeoms',0:tstep:tfmal,yyyO, options); 

out=[yyy(:,2)  yyy(:,3)  yyy(:,4)  yyy(:,5)]; 
fprintf(fde,'%2.8e  \t  %2.8e  \t  %2.8e  \t  %2.8e  \n',ouf); 

%  fprintf(fde,’%2.8e  \t  %2.8e  \t  %2.8e  \t  %2.8e  \t  %2.8e  \n',yyy'); 
fclose(fde); 
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